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About 1% of giants” have anomalously high Li 
abundances in their atmospheres, conflicting directly 
with the prediction of the standard stellar evolution 
models”. This finding makes the production and evo- 
lution of Li in the Universe intriguing, not only in 
the sense of Big Bang nucleosynthesis or the inter- 
stellar medium” , but also for the evolution of stars. 
Decades of efforts have been put into explaining why 
such extreme objects exist’ , yet the origins of Li- 
rich giants are still being debated. Here we report 
the discovery of the most Li-rich giant known to date , 
with a very high Li abundance of 4.51. This rare phe- 
nomenon was observed coincidentally with another 
short-term event: the star is experiencing its lumi- 
nosity bump on the red giant branch. Such a high Li 
abundance indicates that the star might be at the very 
beginning of its Li-rich phase, which provides a great 
opportunity to investigate the origin and evolution of 
Li in the Galaxy. A detailed nuclear simulation is pre- 
sented with up-to-date reaction rates to recreate the 
Li enrichment process in this star. Our results provide 
tight constraints on both observational and theoreti- 
cal points of view, suggesting that low-mass giants can 
internally produce Li to a very high level through 7Be 
transportation during the red giant phase. 


Lithium is too fragile to survive in deeper layers of a 
stellar atmosphere due to the high temperature. Thus 
the first dredge up (FDU) process can sharply dilute the 
surface Li abundance in red giants. This explains why 
the first discovery” of a Li-rich giant evoked great in- 
terests on exploring and understanding the Li-rich ob- 
jects. However, only about 150 Li-rich giants have been 
found PPIE in the past three decades, and ~ 20 
of them were found to be super Li-rich with Li abun- 
dances higher than 3.3. Considering the NLTE correc- 
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tions, three stars were found to be at a level of 


A(Li) > 4.0. Such rare objects could provide a great 
opportunity to reveal the nature of the phenomenon of 
Li-richness because high Li abundance cannot be main- 
tained for a long time due to frequent convection activity. 
Taking advantage of the powerful ability for spectral col- 
lection with the Large Sky Area Multi-Object Fiber Spec- 
troscopy Telescope (LAMOST), we have obtained a large 
sample of Li-rich candidates by measuring the equivalent 
width of the Li ı line at A = 6707.8 Å. One of our candi- 
dates, TYC 429-2097-1, has a super strong Li absorption 
line (see Fig. 1, panel a). We then made a follow-up high- 
resolution observation with the 2.4-m Automated Planet 
Finder Telescope (APF) located at Lick Observatory on 
June 23, 2015. The spectrum covers a wavelength range 
of 374nm — 970nm with a resolution of ~ 80,000. The 
total integration time was 1.5 hours and was divided into 
three single exposures (30 minutes each) for a better sub- 
traction of cosmic-rays. The spectrum of TYC 429-2097-1 
obtained from APF is presented in Fig. 1, panels (b) and 
(e), where the spectrum of HD 48381 is also plotted with 
a vertical shift of +0.3 as a comparison. HD 48381 is a 
star selected from the Gaia-ESO survey DR2, which has 
very similar stellar parameters to TYC 429-2097-1. 


We used the spectroscopic method to derive the stel- 
lar parameters (see the ‘methods’ section for details). 
We presented the final derived parameters of TYC 429- 
2097-1 and the estimated errors in Table 1. The NLTE 
Li abundances for 6707.8 À, 6103.6A, and 8126.3A are 
4.42+0.09, 4.51+0.09, and 4.60+0.08, respectively. The 
averaged Li abundance is A(Li)yrre = 4.51 + 0.09. Com- 
pared to previous studies, TYC 429-2097-1 has the high- 
est Li abundance among all Li-rich giants ever discovered 
(see Fig. 2). The Liabundance in TYC 429-2097-1 is about 
1,000 times as high as the widely-used Li-rich ‘standard’ 
of A(Li)=1.5 (the lower purple dashed line in Fig. 2), 
despite this ‘standard’ being suggested to be luminosity- 
dependent”. It is also about 15 times as high as me- 
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Figure 1: The observed spectra and line profile fittings for TYC 429-2097-1. The low-resolution spectrum from LAMOST is shown in panel 
(a), where Ha line and Lil resonance line are indicated. Panels (b) and (e) display the high-resolution spectrum (blue dot) observed 
by APF near 6708A and 6103 A, respectively. The spectrum of HD 48381 is also plotted (green dots) for comparison in both panel. 
The best profile fittings for these two lines are indicated with the red solid line in panels (c) and (d), respectively, where several 
profiles computed with different Li abundances are also presented. The profile from LTE calculation reaches saturation (cyan curve) 
in panel (c), and increasing Li abundance will barely affect the computed profile, therefore the LTE abundance shown in panel (c) 


is estimated from the equivalent width of this line. 


teoritic Li abundance (the upper purple dashed line in 
Fig. 2), which is thought to be the initial Li abundance 
for newly-formed young stars. 

Although Li-rich giants were reported at various stages, 
such as RGB and core He-burning phases”, the Li-rich 
phase is likely to be a short-term event. An extremely 
Li-rich giant (possibly newly enriched) with rigorous in- 
vestigation on its evolutionary stage would be definitely 
important. The location of the star was derived by the 
maximun likelihood method using the observed parame- 
ters (in this case, Tar, log g, and [Fe/H] derived from the 
spectroscopic method) and a grid of evolutionary mod- 


els computed with the MESA code (see the ‘method?’ sec- 
tion for details). The derived luminosity and mass are 
log(L/L.)=1.95 and M=1.43 Mg, respectively. We used 
the parallax of Gaia DR1™ to test the reliability of the in- 
formation derived from the maximum likelihood method 
independently. The luminosity obtained from Gaia data 
leads to a very similar result of log(Lgai,/Lo)=2.00. The 
mass was tested in the sense that if the mass is well deter- 
mined, the surface gravity from Gaia parallax will show 
good consistency with the spectroscopic log g of 2.25. As 
expected, the final result is log gg,i,= 2.23. Thus, we 
consider that the results derived from the maximum like- 


Table 1: The key information of TYC 429-2097-1 


Property Value 
Name TYC 429-2097-1 
Tote (K) 4696 + 80 
logg 2.25+0.10 
[Fe/H] —0.36 + 0.06 
č: (km s~!) 2.30 + 0.10 
Alli gee 4.51+0.09 
Gaia parallax = (milli-arcsec) 0.73 0.24 
Mass (Mo) 1.434223 
log(L/Lo) 1.95 1529 
log(L¢aia/Lo) 2.00 + 0.06 
log ZGaia 2.23 + 0.16 
bas Oy da © 12.0+ 3.0 
vsini (km s71) 11.341.5 
[a/Fe] 0.19+0.04 


lihood method are reliable, allowing us to robustly lo- 
cate this star on the Hertzsprung-Russel diagram (H-R di- 
agram) along with the corresponding MESA tracks (see 
Supplementary Figure 1). The star is likely occupying the 
region of the RGB-bump, a stage in which the u-barrier 
is destroyed and the enhanced ‘extra-mixing’ might be 
ongoing inside the star. In addition, we also estimated 
the '*C/!SC ratio as it has been suggested that the extra 
mixing will cause a decrease of 1*C/'%C to the range of 
10 — 20. We found that the *C/'°C ratio in this star is 
12.0 + 3.0, which is well within the predicted range. All 
the results obtained above are shown in Table 1. 

It has long been suggested that the Li enrichment 
could be due to contaminations by external sources in 
the enon such as the engulfment of a substel- 
lar component” (e.g. giant planets or brown dwarfs) 
and the accretion from a Li-rich companion or diffuse 
medium. Yet the contribution from external sources is not 
infinite, since the contributor itself has a limited amount 
of Li, typically not higher than 3.3. A simulation on en- 
gulfment of a Jovian planet suggested that a typical upper 
limit for enrichment by such way is ~ 2.27 . Our star has 
a much higher Li abundance than any of those values, 
thus it is very unlikely that the overabundant Li comes 
from the direct contribution of external sources. 

The internal production of Li, on the oair hand, is 
based on the Cameron—Fowler mechanism” (CF mecha- 
nism). The production of ’Be takes place where the tem- 
perature is too high to preserve the newly synthesized ’Li, 
hence 7Be must be transported quickly to the cooler re- 
gion to form Li. This scenario would potentially require 
the low-mass giants to evolve to the RGB-bump, where 
the mean molecular weight discontinuity (or y-barrier, a 
mass gradient caused by FDU) is erased. Meanwhile, it 
would need the presence of deep, enhanced ‘extra mix- 
ing’ (EM) to increase the depth and efficiency of the con- 
vective circulation, which in-turn alters the 12C/'°C into 
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Figure 2: The distribution of Li-rich giants in Tag - Li/H 
plane (left panel) and Tog - A(Li) plane (right 
panel). TYC429-2097-1 is indicated with a red 
star. The cited data are from different litera- 
= as oe in the le if panel, name 

onaco2011 Kirby2012 MS2013 

Adamów2014™, Casey 2016, Balach2000 
RL2005?, Kumar2011™, dd1995™, xr2009™. 
Ruchti2011”, and E All the Li abun- 
dances adopted here are based on NLTE calculations. If 
the original work did not perform the NLTE abundance 
analysis, we appli ied the NLTE corrections interpolated 
from Lind’s grid” to the original LTE abundances. The 
horizontal dashed lines in purple indicate A(Li) =1.5 
and A(Li) =3.3, respectively. The dotted blue line shows 
the upper limit of Li enriched by engulfment of a giant 
planet 


a lower level than that after FDU. The observational fea- 
tures on both the evolutionary stage and *C/C ratio 
of our star coincide with these predictions remarkably 
well, but the limitation of self-production still remains un- 
known in the sense that none of the quantitative calcula- 
tions with a nuclear reaction network has been presented 
to obtain such high amount of Li before. To test this spec- 
ulation, we have made such simulation with a series of pa- 
rameters. By using the RGB stellar structure as the input 
for the EM calculation, with the updated nuclear reaction 
rates and the asymmetric parameters of the EM model, we 
found that A(Li) in the envelope can exceed 4.0 for the 
processed material when the mass circulation finished. 
Our EM calculation with parameters of M = 52, A=0.15, 
fa=0.9 and f,=0.1 (see the ‘methods’ section for the de- 
tails) yields A(Li)= 4.506, where M is the rate of mass 
transport in units of 10° Mg yr, A is log Ty — log T, 

where Ty is the temperature at which the energy released 
from the H-burning shell reaches maximum and T, is the 
maximum temperature sampled by the circulating mate- 
rial, and f4 and f, are the fractional areas of the ‘pipes’ 
occupied by the mass flows moving downward and up- 
ward, respectively, and their values satisfy f4 + f ,=1. This 
reproduces the observed Li abundance for TYC 429-2097- 
lwell. Repeating the same calculation with the alterna- 
tive set of nuclear reaction rates from the JINA database 
yields a similar abundance of A(Li)= 4.515. As a con- 
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Figure 3: The calculated surface abundances and mass fraction of 
3He, ’Be and 7Li as functions of the processing time for 
the mass circulation. The left-hand y-axis indicates the 
mass fractions in logarithmic scale and the right-hand 
y-axis indicates the logarithmic abundances. The initial 
time is set to be zero when the sample material at the 
base of the envelope starts the mass circulation. The 
vertical solid line indicates the boundary between the 
downward and upward motion of the processed mate- 
rial. The yield of "Be reaches maximum by the end of 
the downward motion at t=423 yrs, and then ’Li begins 
to dominate during the upward motion of the sample 
material for 47 yrs. 


trast, assuming this star had never experienced any EM, 
the Li abundance would be constant at the initial value 
of A(Li)=1.16, because the temperature in the envelope 
is too low to ignite both the production and destruction 
reactions of ’Li. The abundances of *He, ’Li and 7Be as 
functions of the processing time for the mass circulation 
are shown in Fig. 3. 

During the EM process, *He is converted to 7Be via 
the reaction of "He(*He, y)’Be, and then 7Be is quickly 
converted to ’Li via the reaction of ’Be(e~, v)’Li. To 
achieve such high level of Li abundance, abundant *He 
is required. The initial surface 3He is computed from the 
MESA model, which is YCHe) = 4.038 x 10-4. Fig. 3 
shows the decrease of *He as a function of the time for EM 
processing. A total amount of Y@He)/H ~ 1.477 x 1074 
is burned-off during this circulation, and the produced 
Y(Li)/H is 3.206 x 1078. This is because another re- 
action, *>He(?He, 2p)*He, dominates over the reaction 
3He(*He, y)’Be, thus consuming the majority of He. The 
strong competition from *He(*He, 2p)*He reaction pre- 
vents more ?He from converting to 7Li. Testing with dif- 
ferent sets of EM parameters shows that the maximum 
of Li abundance from our network calculation is 5.07. 
The *He supply may eventually run out and cannot be 
renewed by the giant, then the surface Li abundance is 
likely to decrease, even if the internal conditions remain 
the same. On the other hand, if the internal conditions do 
change, the surface Li abundance may also decline due 
to the destruction by convective activities in stars. Ei- 
ther way, the super Li-rich phase may disappear after a 


short period of time. In our calculation, the asymmetric 
mass circulation described by a large ratio of f/f, isa key 
factor for achieving super high Li-enrichment. This large 
fa/f, ratio indicates that the upward flow is moving much 
faster (since its ‘pipe’ is thinner) than the downward flow, 
while the mass is conservative in the EM process. 


The cause of the EM has not been well understood, and 
rotationally induced mixing is often attempted. Indeed, 
TYC 429-2097-1 is a slightly rapid rotator with a projec- 
tive velocity of 11.3 kms, which is about ten times faster 
than that for normal giants. The spinning up of an RGB 
star is either caused by the tidal synchronization effects 
in a close binary system or the engulfment of a massive 
planet. We calculated the radial velocities based on 
the two independent observations through LAMOST and 
APF (with an interval of 10 months), and found no sig- 
nificant radial velocity change at a level of a few kilome- 
ters per second, which is the typical uncertainty for the 
RV derived from LAMOST spectra. Thus it is very un- 
likely that a star has a stellar companion which is mas- 
sive and close enough to spin up via tidal synchroniza- 
tion. On the other hand, one would expect some associ- 
ated features that are detectable if a massive planet was 
engulfed and digested. For example, it was found that 
there might be a large probability for Li-rich giants ex- 
hibiting excess in the infrared (IR) flux, yet we found no 
sign of IR excess (see Supplementary Figure 2). In addi- 
tion, if the matter exchange did happen at a certain time, 
there should be some fluctuations in the abundance pat- 
tern. However, TYC 429-2097-1’s '2C/'8C is at a typical 
level for its stage, and its a-abundance is also quite nor- 
mal among the giants with similar [Fe/H]. Given all these 
facts, we speculate that in our case, the enhanced extra 
mixing might neither be caused by the presence of a mas- 
sive planet (if there were any) nor a close stellar compan- 
ion. There are other assumptions often approximated as 
the internal cause of enhanced extra mixing, i.e., thermo- 
haline instabilities and magnetic buoyancy. The thermo- 
haline convection driven by the 7He(#He,2p)*He reaction 
which pene a local depression in the mean molecular 
weight” can cause enhanced extra mixing inside the star. 
The magnetic buoyancy mechanism in the presence of a 
magnetic dynamo would permit the buoyancy of magne- 
tized material near the H-burning shell, thus inducing the 
form of matter circulation in RGB stars”. We speculate 
that the magnetic buoyancy and thermohaline instabili- 
ties might play roles together during the mass circulation, 
in which the former may lead to very fast upward circula- 
tion and the latter drive downward circulation at a much 
slower speed. 


Although the Li abundance measured in this star is su- 
per high, it is still well within the upper-limit that the the- 
oretical model could reach. It is also important to note 
that the RGB-bump is not the only stage for inhabitation 
of Li-rich giants, many Li-rich giants have been reported 
in various stages by previous work, including the core He- 
burning phase, which is very close to the RGB-bump re- 
gion on the H-R diagram. Although it is not preferred by 
our data, if our star occupies this stage, a new scenario 


will be in urgent need for interpreting such high Li abun- 
dance. 
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Methods 


Data Reduction. We followed the standard procedure for data 
reduction with an Interactive Data Language (IDL) package, 
which was originally designed for the FOCES spectrograph : 
The instrumental response and background scatter light were 
considered during the reduction, and cosmic rays and bad pixels 
were removed carefully. The resulting spectrum has a signal-to- 
noise ratio of ~ 160 at 6707.8 Å. 


Deriving the Stellar Parameters. We first combined three iron 
(Fe) line lists PE and calibrated 213 lines out of 257 with the 
solar spectrum . Then we eliminated those seriously blended 
or poorly recognized lines seen from the spectrum of TYC 429- 
2097-1, as well as the lines that are too strong (> 120 mA) or 
too weak (< 20 mA). Finally, 57 Fe 1 and 12 Fe 1i lines are used as 
the parameter indicators. The effective temperature (Tg) is de- 
rived from the excitation equilibrium of Fe 1 lines with excitation 
energy (Eexc) greater than 2.0 v”, The surface gravity (log g) 
is approached by equalizing the two sets of Fe abundances ob- 
tained from Fe 1 and Fe 1 lines, respectively. Statistically, the 
Fe abundance derived from each individual Fe 1 line and the 
equivalent width (EW) from the same Fe 1 line will be mutually 
independent if the micro-turbulence velocity (&,) is correctly 
set. Using this trick, we can obtain €,, and then the metallic- 
ity ([Fe/H]) can be settled simultaneously if all the mentioned 
constraints are achieved. All the Fe abundances are derived from 
NLTE analysis with the MARCS atmospheric models since it 
has been suggested that Fe I lines suffer a non-negligible NLTE 
effect, The procedure of this approach is much more like an 
iteration. We started with the results from LAMOST pipeline as 
the initial input, and then by calculating MARCS models and 
adjusting the stellar parameters step by step, we finally end up 
with a self-consistent solution. Supplementary Figure 3 shows 
the derived Fe abundances from individual lines as functions of 
their EWs (upper panel) and Esx: (lower panel). Based on the 
experience of our previous work using the similar spectroscopic 
method, the errors for Tar, log g, [Fe/H] and €, are estimated to 
be +80K, +0.10 dex, +0.06 dex and +0.10 kms", respectively. 


Determination of the Elemental Abundances. For all the 
species discussed in this paper, we use the spectrum synthe- 
size method to derive their abundances. The theoretical profiles 
of the corresponding lines are calculated based on the MARCS 
model”. An interactive IDL code Spectrum Investigation Util- 
ity (SIU) was applied to calculate the synthetic line profiles. 
The coupled radiative transfer and statistical equilibrium equa- 
tions for the NLTE calculation were solved following the efficient 
method with a revised DETAIL program based on the accelerated 
lambda iteration, we refer readers to Mashonkina et al. (2011) 
for a more detailed description of this method” . The resulting 
departure files are transferred into SIU for NLTE line synthesize. 
The solar iron abundance of log ¢;, = 7.5 was assumed in our 
work. 


In the abundance analysis of Li, the resonance line at 
6707.8A, the subordinate line at 6103.6A and the line at 
8126.3 Aa were used to derive the Li abundance. Although the 
line at 8126.3 A is blended with two telluric lines, it shows sim- 
ilar result to those derived from the resonance and subordinate 
lines. The final Li abundance is determined by averaging the re- 
sults from these three lines. It is noted by many previous studies 
that the NLTE corrections are important for strong lines. In gen- 
eral, the NLTE correction for Li is not large for the ‘Li-normal’ 
stars, however it will significantly increase for Li-rich objects, 


especially for the strong resonance line at 6707.8 A. In the very 
extreme cases (such as ours), the local thermodynamic equilib- 
rium (LTE) theoretical profile of 6707.8A could be saturated at 
the core. Therefore, the NLTE effects were taken into consider- 
ation in our abundance analysis for Li. For the NLTE analysis, 
we applied the same atomic model and line data as those pre- 
sented in Shi et al. 2007. The carbon abundances were de- 
rived from the C1 line at 5086 Å, and the C, line at 5135 a, 
The nitrogen lines were either blended or too weak in our spec- 
trum, so we turned to the CN band near 8003.5 Å to estimate 
the N abundance by fixing C to the value we just derived. Then 
the carbon isotopic ratio was determined by adjusting the con- 
tributions from 1°C and 1°C until we get the best fit to the CN 
band observed profile. The determination of a-abundance (Mg, 
Si, and Ca) with NLTE analysis was based on a series of previous 
work T, The final a-abundance was obtained by averaging 
the abundances obtained from those elements. We also derived 
abundances of several other elements, and they can be found in 
Supplementary Table 1. 

The error of the Li abundance was estimated by changing the 
stellar parameters (namely T.,, log g, and [Fe/H]) within their 
error ranges and calculating the corresponding variations on the 
abundance. The result of this test is presented in Supplementary 
Table 2. It is clear that the Li abundance is more sensitive to the 
variation of T, than that of log g or [Fe/H]. A change of 80 K on 
Teg Will result in a variation of ~ 0.09 dex for the Li abundance. 
Thus, we adopted the variation caused by the error of the ef- 
fective temperature as the uncertainty for each Li I line, which 
is +0.09, +0.09, and +0.08 for the lines of 6707.8 A, 6103.6A, 
and 8126.3A respectively. The error of the final Li abundance 
is obtained by calculating the standard deviation of abundances 
derived from the three lines. For the other elements, if more 
than three lines are used for the abundance determination, we 
calculate their standard deviation and compare it to the error 
caused by the uncertainties of the stellar parameters, and the 
larger one was adopted as the final error. A few of species such 
as *C/+8C ratio were not suitable for the above analysis, we es- 
timated their errors by giving the upper- and lower- limit of the 
best-fit to the profile. 


Maximum Likelihood Method and Evolutionary Stage. The 


likelihood function is expressed following Basu et al.~ , which 
is defined as: 


L= Loosi Lopes Looss us () 
where P »s is the observed parameter (e.g., Teg) and 
1 —(Pobs — Pmodet)” 
Lp, = exp ( obs ; model) : (2) 
j 270 pys 20>... 


The normalized probability of each model p; is expressed as 


4 (3) 
ys L i 
where L; is the likelihood function of the ith model and N,, is 
the total amount of the models. The probability is integrated 
from the boundary constrained by 30 error range of the ob- 
served parameters. Thus, the maximum value of the integrated 
probability is 0.5, and the best-fitted parameters are obtained 
from this probability 

The grid of An cae for calculating the likelihood 
are generated from MESA™ code. The grid covers a wide-range 
of mass from 0.6-3.0 Mẹ with a 0.02 M, interval on mass and a 
0.005 interval on metallicity (Z). The evolution tracks are con- 
structed from the pre-main sequence to the asymptotic giant 


Pi= 


branch (AGB) phase. For generating the grid, the initial pa- 
rameter setup was mostly as same as described in Paxton et al. 
2011 except for the solar chemical abundance (Z/X),. We 
adopted (Z/X), = 0.02294 because calibrated solar model 
with this mixture fits the internal structures from helioseismic 
inversion slightly better than the others” . The MESA p-T 
tables are based on an updated version of Rogers & Nayfonov 
: 

tables™ . Ferguson et al.™ extended the opacity for the so- 
lar composition to the low-temperature case in 2015, and we 
adopted their results in our computation. The stellar metallicity 
was transferred into ‘metal’ abundance Z,,;,, from where the hy- 
drogen and helium abundances (X;,;, and Yjni,) were calculated. 

To obtain the luminosity and log g,,;,, from the parallax, we 
first calculated the bolometric magnitude for the absolute V 
magnitude by Myo) = Ving t BC +5 log 7+5—Ay (in this equation, 
r is used in the unit of arcsec), where the bolometric correction 

i , 

BC was computed following Alonso et al.™ and Ay was esti- 
mated by the Galactic extinction map presented by Schlafly & 
Finkbeiner at 201 12 (all values that were need for the calcula- 
tion were presented in Supplementary Table). We calculated the 


Lo. 

luminosity using the relation of Myo1 — Mbolo = —2.5 log( a ). 
© 

Finally, the log g,,;,was determined following the fundamental 


M T, 
relation logg = log go + log(—) + 4log( =Ë) + 0.4(Mpo — 
M, Teffo 


Mboo), where log go = 4.44, M = 1.36 Mo, Teo = 5777 K, and 
Moo = 4:74mag. The errors were given via the uncertainty 
transfer formula assuming that the errors were contributed by 
the uncertainty of the Gaia parallax. 


Parametric Calculation of Internal Li-enriching. During the 
RGB stage of a low-mass star, we introduced the ‘extra-mixing’ 
or the ‘circulation process’ after the H-burning shell has erased 
the chemical discontinuity left behind by FDU, signatured by 
the bump of the luminosity function in RGB. We followed the 
parameterization of Nollett et a. to perform the parametric 
calculations using M, T,» fa and f, as free parameters. 


The basic assumptions are described as follows. The EM is 
a process of meridional circulation in the radiative zone of low- 
mass red giant stars, the path of mass flow looks like a ‘conveyor 
belt’. A parcel of material with the initial abundance compo- 
sition YF (0) at the base of the convective envelope circulates 
downward through the radiative zone structure and finally re- 
turns to the envelope with newly processed abundance compo- 
sition Ye . The velocity of the sampled material can be expressed 
as dr/dt = M/f4nr?p, where f is equal to f; (or f,) for the 
downward (or upward) circulation, and the density p is a func- 
tion of radius r, governed by the stellar structure. The stellar 
structure, M and A specify as functions of time, the conditions of 
temperature T and density p where the material pass through. 

The computing code integrates the network of reactions of 
the hydrogen burning chain by following the circulation trajec- 
tory. The nuclear reaction rates adopted are from NACRE 
except for two rates, the updated rate of ’Be(p,y)*B is from 
Du et al. 2015 ~ and the rate of ’Be(e, v)’Li from the JINA 
database. The convective overturn-time of the envelope is set to 
~1 yr, which is a good approximation in the sense that the mix- 
ing between the processed material and the convective envelope 
is instantaneous. The abundance of the i-th nucleus changes in 
the envelope during the transport is due to the nucleosynthe- 
sis and the material replacement, and this corresponds to inte- 
grating Y,° = (M/M; )(Y? — YE), where M; is the mass of 
convective envelope. In order to obtain the RGB stellar struc- 
ture, including in particular the initial abundance composition 


in the envelope and the distributions of temperature and den- 
sity in the radiative zone, we calculated the stellar evolution 
model for M = 1.36 M, by the MESA code. The initial abun- 
dance of ’Li in the sample material at the base of the envelope 
is ~ 1.024 x 1071 (A(Li)=1.86) obtained from the present stel- 
lar evolutionary model calculation, and the 7Li abundance may 
increase to a level of A(Li) exceeding 4 in the processed material 
when the mass circulation finishes. 

During the downward mass circulation, the abundance 
of 7Be increases quickly because the construction reac- 
tion *He(*He,7)’Be wins against the destruction reaction 
7Be(p, y)°B, and the maximum yield of ’Be is around the turn- 
ing point of the mass circulation, where the temperature reaches 
the highest value T,. In contrast, the production change of 7Li 
behaves rather dramatically during the downward mass circu- 
lation due to the complex competition between the production 
reaction ’Be(e~, v)’Li and the destruction reactions ’Li(p, y)*®Be 
and ’Li(p, a)*He. The abundance of ’Li drops suddenly at about 
200 yrs of the processing time as the rates of destruction reac- 
tions increase quickly with increasing temperature, and the 7Li 
abundance keeps very low during the downward mass circula- 
tion although it increases slightly after about 230 yrs due to the 
decay of fast growing abundant Be. In contrast, the abundance 
of ’Li increases sharply during the upward mass circulation due 
to the fast decreasing of the destruction reaction rates of 7Li 
as the temperature decreases quickly. The processed abundance 
Y/ finally reaches a super-high-level of A(Li)= 4.506. The abun- 
dance of ’Li in the envelope contains the contributions of the 
mass circulation and mass replacement processes, and the to- 
tal processing time can be estimated as M,/M, which is about 
2.1 x 10t yrs by using the average value of the envelope mass 
for this RGB star. 


The Projected Rotational Velocity. Following the assumption 


of Bruntt et al. the external broadening of the line profile 
was assumed to be contributed from the stellar rotation, 
instrumental broadening, and macro-turbulence. The projected 
rotational velocity (vsini) was derived by using five isolated 
iron lines at 6151A, 6229A, 6380A, 6703A, and 6810A. 
The instrumental broadening was calculated from fitting the 
emission lines of the arc lamp with a Gaussian profile. The 
macro-turbulence velocity was estimated using the relation of 
Hekker et al. which is a function of T. and logg. Then 
we calculated a set of the theoretical spectra broadened with 
different rotational velocities, and vsini was determined by 
finding the best-fit to the observed iron line profiles. 
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Supplementary Information 


Supplementary Table 1. Other information of TYC 429-2097-1 


Property Value 
Position (J2000) R.A. 17:53 : 46.07 
DEC. 06 : 42 : 41.20 
l 32.4111 
b 15.8647 
Vinge (mag) 11.27 (UCAC4) 
Gaia parallax = (milli-arcsec) 0.73 0.24 
Ay (mag) 0.47 
BC (mag) —0.385 
Mhol (mag) —0.27 
[C/Fe] —0.02 + 0.10 
[N/Fe] 0.45 0.10 
[Mg/Fe] 0.23 + 0.05 
[Si/Fe] 0.19 + 0.01 
[Ca/Fe] 0.16 + 0.06 


Supplementary Table 2. The variation of Li abundance caused 
by the uncertainties of stellar parame- 


ters 
Wavelength AT. Alogg  A[Fe/H] 
80K 0.10dex 0.06dex 
6103.6A 0.09 0.02 0.07 
6707.8A 0.09 0.03 0.06 
8126.3A 0.08 0.01 0.07 
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Supplementary Figure 1. The corresponding evolutionary 


tracks computed from the MESA 
code for TYC 429-2097-1. The mass 
range of the tracks is from 1.0 M, 
(bottom) to 1.8M, (top), with an 
interval of 0.2 Mọ. The evolutionary 
stages for each track are indicated 
with different colors, namely: ‘black’ 
— SGB, ‘dark blue’ — RGB, ‘red’ — 
RGB-bump and ‘orange’ — RGB tip. 
The position of TYC429-2097-1 is 
indicated with a black dot. The error 
bars indicate the uncertainties in Tar 
and logg, which are adopted to be 
+80K and 0.10 dex from our stellar 
parameters determination method. 
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Supplementary Figure 2. 


The IR-excess diagram for TYC 429- 
2097-1. TYC429-2097-1 is is indi- 
cated with a red dot, while the other 
stars are from Rebull et al. 2015. 
The horizontal line at A(Li) = 1.5 
shows the adopted division between 
Li-rich and normal stars, and the ver- 
tical line at [3.4] — [22] = 0 indi- 
cates the photospheric locus. The blue 
dots denote the objects with strong IR- 
excess, while the black dots represent 
the stars with no significant IR-excess. 
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Supplementary Figure 3. The determination of the stellar pa- 


rameters. The figure shows the ab- 
solute NLTE abundance from lines of 
Fe 1 (red open circle) and Fe 1 (blue 
dots) in TYC 429-2097-1 as functions 
of their EWs (top panel) and Exe 
(bottom panel). The slopes are indi- 
cated in the corresponding panel. The 
mean A(Fe) averaged from Fe 1 and 
Fe 11 lines is shown with a dash-dotted 
line in both panels. The final stel- 
lar parameters are denoted in the se- 
gence of Tag (K), log g, [Fe/H], and €, 
(kms~') in the bottom panel. 
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